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COORDINATE-INVARIANT INCREMENTAL LYAPUNOV FUNCTIONS 

MAJID ZAMANI^ AND RUPAK MAJUMDAR^ 

Abstract. The notion of incremental stability was proposed by several researchers as a strong property of 

dynamical and control systems. In this type of stability, trajectories converge to each other, rather than to an 

equilibrium point or a particular trajectory. Similarly to stability, Lyapunov functions play an important role 

^ ' in the study of incremental stability. In this paper, we propose coordinate-invariant notions of incremental 

Q ' Lyapunov function and provide characterizations of incremental stability in terms of existence of the proposed 

\—^ ^ Lyapunov functions. 

in 

u 

. , Incremental stability requires that all trajectories converge to each other, rather than to an equilibrium point 

i^H ' or a particular trajectory. While it is well-known that for linear systems incremental stability is equivalent to 

C^ , stability |Ang02| , it can be a much stronger property than stability for nonlinear systems. See [ZTllj for a 

historical discussion and a broad list of applications of incremental stability. 

Similarly to stability, Lyapunov functions play an important role in the study of incremental stability. In 
^S| ■ [Ang02| , Angeli proposed the notions of incremental Lyapunov function and incremental input-to-state Lya- 

^ I punov function, and used these notions to prove characterizations of incremental global asymptotic stability 

"J ■ ((5-GAS) and incremental input-to-state stability (^-ISS). Notions of (5-GAS, (5-ISS and incremental Lyapunov 

1 (»«. I functions, proposed in [Ang02, , are not coordinate invariant, in general. Since most of the controller design 

^vq ■ approaches benefit from changes of coordinates, in |ZT11| . authors proposed different notions of J- GAS and 

<5-ISS which are coordinate invariant. In this paper, we propose notions of incremental Lyapunov function 
i— N . and incremental input-to-state Lyapunov function that are coordinate invariant. Moreover, we use these 

new notions of Lyapunov functions to characterize the notions of incremental stability, proposed in |ZT11| . 

In [ZMllj . we have used coordinate invariant incremental Lyapunov functions as a basis for the design of 

controllers enforcing control systems incrementally stable. 

H 2. Control Systems and Stability Notions 

2.1. Notation. The symbols N, M, M''' and M(j" denote the set of natural, real, positive, and nonnegative 
real numbers, respectively. Given a vector x G K", we denote by Xi the i-th element of a;, and by ||a;|| the 
Euclidean norm of x; we recall that ||a;|| = \J x\ + x"^ + ■■■ + x"^. Given a measurable function / : Rq — > R", 
the (essential) supremum of/ is denoted by l|/||oo; we recall that ||/||oo '■— (ess)sup{||/(i)||,i > 0}. Function / 
is essentially bounded if ||/||oo < oo. For a given time r G M+, define fr so that frit) — f{t), for any t G [0, r), 
and f{t) = elsewhere; / is said to be locally essentially bounded if for any r e M+, fr is essentially bounded. 
A continuous function 7 : Rq -^ Rj, is said to belong to class JC if it is strictly increasing and 7(0) = 0; 7 is 
said to belong to class /Coo if 7 G /C and 7(r) — > 00 as r — > 00. A continuous function /3 : Rq x Rq — >■ Rq is 
said to belong to class ICC if, for each fixed s, the map /?(r, s) belongs to class /Coo with respect to r and, for 
each fixed nonzero r, the map /3(r, s) is decreasing with respect to s and /3(r, s) — > as s — ^ 00. A function 
d : K" X R" —i- R'^ is a metric on E" if for any x,y,z G M", the following three conditions are satisfied: i) 
d{x,y) = if and only ii x — y, ii) d{x,y) — d{y,x); and iii) d{x, z) < d{x,y) + d{y,z). For a set ^ C M", 
and any x G M", d(a;,^) denotes the point-to-set distance, defined by d(a;,^) — miy^j,d{x,y). 
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2.2. Control Systems. The class of control systems that we consider in this paper is formalized in the 
following definition. 

Definition 2.1. A control system is a quadruple: 

S-(M",U,W,/), 
where: 

• M" is the state space; 

• U C M™ is the input set; 

• U is the set of all measurable functions of time from intervals of the form ]a, 6[C M to U with a < 
and b > 0; 

• / : R" X U — > M" is a continuous map satisfying the following Lipschitz assumption: for every compact 
set Q C M", there exists a constant Z 6 K+ such that \\f{x,u) — /(?;,u)|| < Z\\x — y\\ for all x,y (£ Q 
and all u E U. 

A curve £, :]a, 6[-> M" is said to be a trajectory of E if there exists v €U satisfying: 

(2.1) m^fim.Ht)), 

for almost all t £ ]a,b[. We also write £,xv(t) to denote the point reached at time t under the input v from 
initial condition x — S,xviO)', this point is uniquely determined, since the assumptions on / ensure existence 
and uniqueness of trajectories jSon98j . A control system S is said to be forward complete if every trajectory 
is defined on an interval of the form ]a, oo[. Sufficient and necessary conditions for a system to be forward 
complete can be found in |AS99j . A control system E is said to be smooth if / is an infinitely differentiable 
function of its arguments. 

2.3. Stability notions. Here, we recall the notions of incremental global asymptotic stability ((J^-GAS) and 
incremental input-to-state stability ((5g-ISS), presented in jZTllj . 

Definition 2.2 ( |ZTllj ). A control system E is incrementally globally asymptotically stable (S^-GAS) if it is 
forward complete and there exist a metric d and a JCC function /3 such that for any t G Rq , any x,x' e R" 
and any v ElA the following condition is satisfied: 

(2.2) d(C,4i),C.'..(t))</5(d(x,x'),t). 

As defined in |Ang02| , (5-GAS requires the metric d to be the Euclidean metric. However, Definition 12.21 onlv 
requires the existence of a metric. We note that while J-GAS is not generally invariant under changes of 
coordinates, Jg-GAS is. 

Definition 2.3 ( [ZTll] ). A control system S is incrementally input-to- state stable (5^-ISS) if it is forward 
complete and there exist a metric d, a ICC function (3, and a /Coo function 7 such that for any t e M.q , any 
x,x' G R", and any v, v' £U the following condition is satisfied: 

(2.3) d {^xvit)A.'v'{t)) < /3 (d {x, x') , i) + 7 (11^ - v'WJ . 



By observing p.2|) and (|2.3p . it is readily seen that (5g-ISS implies (5g-GAS while the converse is not true in 
general. Moreover, whenever the metric d is the Euclidean metric, (5g-ISS becomes S-ISS as defined in |Ang02| . 
We note that while S-ISS is not generally invariant under changes of coordinates, (5g-ISS is. 

Here, we introduce the following definition which was inspired by the notion of uniform global asymptotic 
stability with respect to sets in [LSW96 . 

Definition 2.4. A control system E is uniformly globally asymptotically stable (U^GAS) with respect to a set 
A if it is forward complete and there exist a metric d, and a ICC function (3 such that for any t G Rq , any 
x G R" and any v £ U the following condition is satisfied: 

(2.4) d{^xAt),A)<l3{d{x,A),t). 
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We discuss in the next section characterizations of iJ^-GAS and (J^-ISS in terms of existence of incremental 
Lyapunov functions. 

2.4. Characterizations of incremental stability. This section contains characterizations of (5g-GAS and 
(Jg-ISS in terms of existence of incremental Lyapunov functions. We start by defining the new notions of 
(53-GAS and 5g-ISS Lyapunov functions. 

Definition 2.5. Consider a control system S = (M", \J,iJ, f) and a smooth function V : M" x M" -^ M.q . 
Function V is called a Ss-GAS Lyapunov function for S, if there exist a metric d, /Coo functions a, a, and 
K G IR+ such that: 

(i) for any x, x' G R" 

a{d{x,x')) < V{x,x') <a{d{x,x')); 
(ii) for any x, x' G R" and any m G U 

^f{x,u) + §^fix',u)<~KVix,x'). 

Function V is called a 6j-ISS Lyapunov function for S, if there exist a metric d, IC^o functions a, a, a , and 
K G R^ satisfying conditions (i) and: 



(iii) for any x, x' G R" and for any u, m' G U 

11/(2^, «) + If/(^', "') < -^V{x, x') + a{\\u - u'W). 
Remark 2.6. Condition (iii) of Definition \2.5\ can be replaced by: 

^f{x, u) + |^/(x', u') < -p{d{x, x')) + a{\\u~ u'W), 

where p is a /Coo function. It is known that there is no loss of generality in considering p(d(x,x')) = nV{x,y), 
by appropriately modifying the S^-ISS Lyapunov function V (see Lemma 11 in JPW96J ). 

While (5-GAS and (5-ISS Lyapunov functions, as defined in |Ang02| , require the metric d in condition (i) in 
Definition l2.5l to be the Euclidean metric, Definition 12. 51 only requires the existence of a metric. We note that 
while (5-GAS and (5-ISS Lyapunov functions are not invariant under changes of coordinates in general, Jg-GAS 
and (5g-ISS Lyapunov functions are. 

We now introduce the following definition which was inspired by the notion of uniform global asymptotic 
stability (UGAS) Lyapunov function in |LSW96j . 

Definition 2.7. Consider a control system S, a set A, and a smooth function V : M" — > R(j". Function V is 
called a U^GAS Lyapunov function, with respect to A, for S, if there exist a metric d, /Coo functions a, a, 
and K G R^ such that: 

(i) for any x G R" 

a{d(x,A)) < V{x) <a{d{x,A)); 
(ii) for any x G R" and any u G U 

^f{x,u)<~nV{x). 

The following theorem characterizes UgGAS in terms of existence of a UgGAS Lyapunov function. 
Theorem 2.8. Consider a control system S and a set A. If \J is compact and d is a metric such that the 



SU! 



function ^ji^) — d(a;,y) is continuous^ for any y G R" then the following statements are equivalent: 

(1) S is forward complete and there exists a U^GAS Lyapunov function with respect to A, equipped with 
the metric d. 

(2) S is UjCAS with respect to A, equipped with the metric d. 



Here, continuity is understood with respect to the Euchdean metric. 
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Proof. First we show that the function (t>{x) = d(x,A) is a continuous function with respect to the Euchdean 
metric. Assume {a^nj^i is a converging sequence in M" with respect to the Euchdean metric, implying: 
a;„ — )> X* as n — >■ oo for some x* e M". By triangle inequality, we have: 

(2.5) d{x*,y) <d{x*,xn) + d{y,xn), 
for any n G N and y £ A. Using inequality (|2.5p . we obtain: 

(2.6) cj)ix*)^ Md{x*,y) < inf {d (x*,a:„) + d (2/,x„)} 

yeA yeA 

= inf d(y, a;„) +d(x*,a;„) 
yeA 

= (j) {x„) + d {x* , x„) . 
Using inequality (|2.6p and the continuity assumption on d, we obtain: 

(2.7) (/)(x*) < lim inf<?!)(a;„), 

n— foo 

for any n £ N, where limit inferior exists because of greatest lower bound property of real numbers [RRA09) . 
By doing the same analysis, we have: 

(2.8) (j){x*) > lim sup</)(x„), 

n— >oo 

for any n G N. Using inequalities (|2.7p and (|2.8p . we obtain: 

(2.9) </'(a;*)= lim 0(x„), 

implying that is a continuous function. Since (p^x) = d{x,A) is a continuous function, by choosing 
uJi{x) — uJ2{x) — d{x,A) and using Theorem 1 in |TP00j . the proof completes. D 

Before showing the main results, we need the following technical lemma, inspired by Lemma 2.3 in | Ang02] . 
Lemma 2.9. Consider a control system S = (M",U,W, /). // S is 6^-GAS, then the control system S 

diagonal set A, defined by: 
(2.10) 



U,iJ,f), where fiCjv) = [f {£,1 , v)'^ , f {^2 , v)'^] , and C, = [^'[,^2] > '■' U^GAS with respect to the 



A = {ze 



pxe 



Proof. Since E is (5g-GAS, there exists a metric d : . 
we define a new metric d 

(2.11) 



X 



M.Q such that property (12.21) is satisfied. Now 



-^ K+ by: 

d(z, z') = d{xi,x'^) + d(a;2, 2^2), 

E M^". It can be readily checked that d satisfies all 



,T iT 
Xi , X2 



for any z = [xi'^,X2'^] S M^" and z' = 

three conditions of a metric. Now we need to show that d(z. A), for any z — [xj, x^l G R^", is proportional 
to d(xi,a;2). We have: 



(2.12) 



d{z,A) 



inf d(z, z') = inf d 
z'eA x'eR" 



Xl 
X2 



X 

x' 



= inf (d(a;i,x') +d(x2,2;')) < d(a;i,X2) 



Since d is a metric, by using the triangle inequality, we have: d(a;i, 2:2) < d(a;i, x') + d(x2, x') for any x' G R", 
implying that d(a;i,X2) < d(z. A). Hence, using (|2.12l) . one obtains: 



(2.13) 



d(xi,a;2) < d(z. A) < d(xi,a;2) =^ d(a;i,a:2) = d(z. A). 
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Using equality (|2.13p and property (I2.2p . we have: 

(2.14) d(C...(t),A) - d(C.,„(t),C...(i)) 

< /3(d(a:i,X2),t) = /3(d(z,A),t), 

for any t G M(J", and v eU, where (^v — [^xivjix2v] j ^^d z = [a;f,a;^] . Hence, S is UgGAS with respect 

to A. D 

We can now state one of the main results, providing characterization of (5g-GAS in terms of existence of a 
(Jg-GAS Lyapunov function. 

Theorem 2.10. Consider a control system E. // U is compact and d is a metric such that the function 
ip{x) = d{x,y) is continuouJQ for any y G M" then the following statements are equivalent: 

(1) E is forward complete and there exists a 5j-GAS Lyapunov function, equipped with the metric d. 

(2) E is S^-GAS, equipped with the metric d. 

Proof. The proof from (1) to (2) has been showed in Theorem 2.6 in |ZM11) . even in the absence of the 
compactness and continuity assumptions on U and d, respectively. We now prove that (2) implies (1). Since 
E is Jg-GAS, using Lemma [^31 we conclude that the control system E, defined in Lemma [^7^ is U3GAS with 
respect to the diagonal set A. Since il^{x) — d{x,y) is continuous for any y G R", it can be easily verified 
that ip{z) — d(z, z') is also continuous for any z' G M^", where the metric d was defined in Lemma [2791 Using 
Theorem 12.81 we conclude that there exists a UgGAS Lyapunov function V : M'^" — > Rq , with respect to A, 
for E. Thanks to the special form of E and using the equality (|2.13p . the function V satisfies: 



(i) a{d{x,x')) < V{x,x') <a{d{x,x')); 
(ii) ^f{x,u) + §^f{x',u)<-KV{x,x'), 

for some /Coo functions a, a and some k G M+. Hence, V is a (5g-GAS Lyapunov function for E. D 



Before providing characterization of Jg-ISS in terms of existence of a (5g-ISS Lyapunov function, we need the 
following technical lemma, inspired by Proposition 5.3 in |Ang02| . To state the following results, we need to 
define the function: 

u if w G U, 



(2.15) satu(M) - ^1 . II , II -r ^ 1 1 

\ ' 1 argmmti'gu ||m — w|| liu^V. 

As explained in | Ang02| , by assuming U is closed and convex and since || • || : K™ — > Rq is a proper, convex 
function, the definition (12.15^ is well-defined and the minimizer of \\u' — u\\ with w' G U is unique. Moreover, 
by convexity of U we have: 

(2.16) \\satu{ui)-satu{u2)\\<\\ui~U2l \fui,U2£W. 

Lemma 2.11. Consider a control system E = (R",U,i^, /), where U is closed and convex. If Yi is S^-ISS, 
equipped with a metric d such that ip{x) = d{x,y) is continuous for any y G R", then there exists a /Coo 
function p such that the control system E ~ (R^", D, 2?, /]f| is U^GAS with respect to the diagonal set A, 
where: 

/(Ci,satu(wi+p(d(Ci,6))'^2)) 
/(6,satu(tJi -p(d(Ci,6))'^2)) 



(2.17) /(Cc.) 



C = [^^,^2^^]^, D = U X Z?i(0), andu;^ [u^I^u;^]'^- 



Here, continuity is understood with respect to the Euchdean metic. 

V is the set of aU measurable, locaUy essentially bounded functions of time from intervals of the form ]a, 6[C R to D with 
a < and b > 0. 
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Proof. The proof was inspired by the proof of Proposition 5.3 in |Ang02| . Since S is (J^-ISS, equipped with 
the metric d, there exists some ICC function (3 and /Coo function 7 such that: 

(2.18) d(e.,(t),Cx'^'(i)) <max{/3(d(x,x'),t),7(l|t^-w'lloo)}. 

Note that inequahty (|2.18p is a straightforward consequence of inequahty p.3p (see Remark 2.5 in }SW95) ). 
Using Lemma [2.91 and the proposed metric d in p.ll[) . we have: d(x,x') = d(z, A), where z — [x^,a:'"^] . 
Without loss of generahty we can assume a{r) = (3{r, 0) > r for any r g M+. Let p be a /Coo function satisfying 
p{t) ^ l7~^ ° (q^~"'^('')/4)- Now we show that 



7 (\\2u;2{t)p (d(C..(t), A)) II) < d(z, A)/2 



^0 ' 



any z G K^", and any uj £ V. Since 7 is a /Coo function and uj2{t) G 'Bi(O), it is enough to 



7 2p d(C.^(t),A) <d(z,A)/2 



(2.19) 

for any t G 
show 

(2.20) 
Since 

(2.21) 7(2p(d(C..(0),A))) -7(2p(d(z,A))) < q-^ (d(z, A)) /4 < d(z, A)/4, 

and (p{z) = d(z, A) is a continuous function (see proof of Theorem 12. 8p . then for aU t G M.^ small enough, we 
have 7 (2p (d(C2„(t), A))) < d(z, A)/4. Now, let 

(2.22) h=mi{t>0 \j{2p(d{CUt),^))) >d(z,A)/2}. 

Clearly ii > 0. We will show that ti = 00. Now, assume by contradiction that ti < 00. Therefore, the 
inequality (|2.20p holds for all t G [0,ti). Hence, for almost all t G [0,ii), one obtains: 

(2.23) j(J2uj2{t)p{d{C.Ut),^))\\) <7(2p(d(Gc.(i),A))) < d(z, A)/2 < a (d(z. A)) /2. 
Let V and v' be defined as: 

v{t) = satu (uji{t) + p (d(G^(i), A)) W2(i)) , 

v'{t) = sat^JLi{t)-p(d{CUt),^))^2{t)] 



By using ((2A6| . we obtain: \\v{t) - v'{t)\\ < 2uj2{t)p ('d(G^(i), A)) . Using (|2J8)) and ([2^, we have 



(2.24) d{CUt). A) = d (C,„(t), C.'.' (t)) < /3 (d(x, x'), 0) = /3 (^d(z. A), Oj = a {d{z, A)^ 

for any i G [0,ii] and any z = [a;^,x'^]^ G M^" which implies that 7 [2^ Cd(Gc^(t), A))) < d(z,A)/4, 

contradicting the definition of ti. Therefore, ti = 00 and inequality (|2.19l) is proved for all t G M.^ . Therefore, 
using (PIB and (PTTO)) . we obtain: 

(2.25) d(C.c.(t),A)-d(e.,(t),e.'.,(i)) < max{/3(d(a:,a;'),<),7(lk-t^'lloo)} 

< max {/3 (d(a;, a;'), i) , 7 (||2w2P (d(C.c., A)) ||^) } 

< max|/3('d(z,A),tVd(z,A)/2|, 



for any z = [x , x \ G M ", any oj G 2?, and any i G KJ. Since /? is a fCC function, it can be readily seen 
that for each r > if d(z. A) < r, then there exists some T^ > such that for any t > T^, /3 I d(z. A), i j < r/2 

and, hence, d(C2^(i), A) < r/2. For any e G M+, let fc be a positive integer such that 2^^r < e. Let r\ — r 

and r, = rj_i/2 for i > 2, and let t ^Tri+Tr2^ h T^fc. Then, for t > t, we have d{C^^{t),A) < 2-^r < e 

for all d(z. A) < r, all u G V, and all t > t. Therefore, it can be concluded that the set A is a uniform 
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global attractor for the control system E. Furthermore, smce d{(^z(^{t),A) < j3 (d(z, A),0) for all t e M([, 

all z G M^", and all uj ^ V, the control system S is uniformly globally stable and as showed in [TPOOj . it is 
UgGAS. D 

The next theorem provide characterization of (5g-ISS in terms of existence of a (5g-ISS Lyapunov function. 

Theorem 2.12. Consider a control system E. If \J is compact and convex and d is a metric such that the 
function V'(a^) = d(x,2/) is continuous^ for any y £ R" then the following statements are equivalent: 

(1) E is forward complete and there exists a S^-ISS Lyapunov function, equipped with metric d. 

(2) E is Sj-ISS, equipped with metric d. 

Proof. The proof from (1) to (2) has been showed in Theorem 2.6 in |ZM11) . even in the absence of the 
compactness and convexity assumptions on U and the continuity assumption on d. We now prove that (2) 
implies (1). As we proved in Lemma 12.111 since E is (5g-ISS, it implies that the control system E, defined 
in Lemma [2.111 is UgGAS. Since ip{x) = d{x,y) is continuous for any y € R", it can be easily verified that 
■0(2;) = d{z, z') is continuous for any z' G R^", where the metric d was defined in Lemma [Z!9l Using Theorem 
12. 8| we conclude that there exists a UgGAS Lyapunov function V, with respect to A, for E. Thanks to the 
special form of E and using the equality (|2.13p . the function V satisfies: 

(2.26) a(d(x, x')) < V{x, x') < aid{x, x')), 
for some /Coo functions a, a, any x,x' £ R", and 

dV dV 

(2.27) —fix,sat^){di + p{d{x,x')))d2) + g^fix',ss>tuidi - p{d{x,x'))d2)) < -kV{x,x'), 

for some k e R+ and any [df ,^1^] G D. By choosing di — {u + u')/2 and d2 — {u — u')/{2p{d{x,x'))) for 

any u, u' £ U, it can be readily checked that [d[ , d^] G U x Bi{0), whenever 2p(d(a;, x')) > \\u — u'\\. Hence, 
using p.27p . we have : 

dV dV 

(2.28) V>{d{x,x')) > \\u-u'\\ ^ —f{x,u) + —fix',u') < -kV{x,x'), 

where ip{r) = 2,p(r). As showed in Remark 2.4 in |SW95j . there is no loss of generality in modifying inequality 
((2:281) to 

dV dV 

(2.29) —f{x, u) + —fix', u') < ~W{x, x') + 7(||k - u'\\), 

for some /Cqo function 7 and some k G R^, which completes the proof. D 
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